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' We present a detailed theoretical study on the thermodynamic properties of impure quasi-one 

^vQ , dimensional unconventional charge-, and spin-density waves in the framework of mean- field theory. 

The impurities are of the ordinary non-magnetic type. Making use of the full self-energy that 
takes into account all ladder-, and rainbow- type diagrams, we are able to calculate the relevant low 
temperature quantities for arbitrary scattering rates. These are the density of states, specific heat 
and the shift in the chemical potential. Our results therefore cover the whole parameter space: they 
include both the self-consistent Born and the resonant unitary limits, and most importantly give 
' exact results in between. 
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I. INTRODUCTION 



In the past few decades many papers have been published dealing with the mean-field theoretical study of impurity 
scattering in either density waves (DW) or superconductors (SC). It has been pointed out earlier by several authors^ 2 -^ 
that non-magnetic impurities should have a pair-breaking effect on conventional density waves, similar to the one 
caused by paramagnetic impurities in SCs, found first by Abrikosov and Gor'kov (AG) A In this connection, the 
original AG theory;^ developed specifically for the treatment of dilute, randomly distributed magnetic impurities 
' in s-wave BCS superconductors, was later successfully applied with minor modifications for example in excitonic 
insulators^ and in the context of conventional charge and spin density waves with constant gap3^ All these studies 

■ however, particularly because of the assumed dilute concentration, considered the impurity effects on the level of 
second order self-consistent Born approximation, and retained only the first nontrivial self energy insertion in the 

qq ■ proper self energy caused by impurities. On the other hand, in case of strong scatterers, higher-order diagrams in the 
f**«« ' self energy might be of importance, especially at low temperatures in view of the anomalous behavior of the normal 
lO '. state (Kondo effect) 

After the discovery of unconventional superconductors, the extension of the field of density waves to DWs with 
wave- vector dependent gap A(k) (termed unconventional) arose in a natural way. Unconventional density waves 
(UDW), cither spin or charge, have the interesting property, that though a robust thermodynamic phase transition 
is clearly seen in experiments, they lack any spatial modulation of charge or spin density due to the zero average of 

■ the gap over the Fermi surface. Consequently, the phase transition cannot be caught in the act with conventional 
means like x-ray or NMR. This state of affairs is usually referred to as the "hidden order" in recent literature J£ One 
of the main reasons of interest on UDWs arises from high-T c superconductors, where one of the competing models 
in the pseudogap phase is the e?-density wave state, a kind of two-dimensional unconventional charge density wave 
(UCDW)i 14 ' 15 i 16 The relevance of UDWs in general is further emphasized by recent studies such as: the mysterious 
low temperature micromagnetic phase of some heavy fermion compounds seem to be described rather well by un- 
conventional spin density wave (USDW), 17 -^ the pseudogap phase of (TaSe^I is explained, at least on qualitative 
level, by a coexisting CDW+UCDW order ; 19 ' 20 and finally the ground state of the transition metal dichalcogenid 
2H-TaSe2 was attributed to an /-wave UCDW lately^ As to the quasi-one dimension, which is the natural place 
of occurence of DWs, one of the most likely candidate for possessing some kind of CDW of unconventional type in 
its low temperature phase is the organic conductor a-(BEDT-TTF)2KHg(SCN)4. In the past few years this salt has 
been investigated extensively ; 22 ' 23 ' 24 ' 25 ' 26 ' 27 and recently the experimental findings of these latter works have been 
explained rather well by a quasi-one dimensional UCDW, including the threshold electric field ) 28 ' 29 ! 30 the angular 
dependent magnetoresistance) 31 ' 32 and the magnetothermopower and Nernst effect. 

Impurity scattering in UDWs and unconventional SCs has been investigated many times from different aspects, but 
the theoretical formulations were limited to include either the weak scattering Born or the resonant unitary limit i 34 ' 35 ' 36 
The main outcome of these studies is that the presence of non-magnetic impurities leads to the destruction of ground 
state. The present work is partly motivated by the surprising discovery that the Ni impurities enhance the energy gap 
of d-wave DW in underdoped YBCO. In light of this, our aim in this paper is to explore the intermediate regime 
as well by making use of the full self energy in the non-crossing approximation. We give a detailed analysis of the 
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thermodynamic properties of UDWs, without making any constraint on the strength of the forward and backscattering 
parameters. As we shall see, our general formulas suggest that non-magnetic impurities always suppress density wave 
ground state. The enhancement of the pseudogap in YBCO by magnetic Ni impurities is proposed to originate from a 
different mechanism discussed in a separate publication, 3 - Calculations similar to this has been carried out for heavy 
fermion superconductors 39 and d-density wave s ' . To the best of our knowledge however, this question in quasi-one 
dimensional UDWs has not been addressed so far. 

The article is organized as follows: in Section 2 we present the formalism necessary to treat impurity scattering in the 
non-crossing approximation in unconventional density waves. In Section 3 we calculate the relevant thermodynamic 
quantities in the low temperature phase. Section 4 is dedicated to the analysis of the quasiparticle density of states, 
and finally Section 5 is devoted to our conclusions. 



II. FORMALISM 

We begin with the mean-field Hamiltonian describing pure quasi-one dimensional density wave s 42 ' 43 

#o = £( k )( a t* a ^v - a k-Q. CT a k-Q.a) + A (T (k)a+ (T a k _Q !CT + A* (k)a^_ Q)0 .a k)(7 , (1) 

k,(T 

where aJ CT and ak, CT are the creation and annihilation operators for an electron in a single band with momentum k 
and spin a. The best nesting vector is given by Q = (2fep, 7r/6, 7t/c), and the k sum in Eq. (fl]) is restricted to the 
reduced Brillouin zone as usual: k x runs from to 2kp, with kp being the Fermi wave vector. Our system is based 
on an orthorhombic lattice with lattice constants a, 6, c towards directions x,y,z. The system is highly anisotropic, 
the kinetic energy spectrum of the particles linearized around the Fermi surface reads as 

£(k) = vp(k x — hp) — 2tt,cos(bk y ) — 2t c cos(ck z ), (2) 

and the quasi-one dimensional direction is the x axis. For (U)CDWs the density wave order parameter A CT (k) is even, 
while for (U)SDW systems it is an odd function of the spin index. In the case of conventional density waves the order 
parameter is constant on the Fermi surface^ 2 - As opposed to this, for unconventional condensates it depends on the 
perpendicular momentum and has different values at different points on the Fermi surface. The precise k-dependence 
is determined by the matrix element of the electron-electron interaction through the gap equation^ 3 - 
The interaction of the electrons with non-magnetic impurities is described by the Hamiltonia n 34 ! 35 ' 44 

Hi = ± E e-^^k + qWR,-)*^), (3) 

k,q,<T,j 

U( ^>-{u*(Q)e^ 17(0) J' (4) 

where V is the sample volume, Rj is the position of the jth impurity atom on the lattice, and the momentum space 
is divided into k and k — Q subspaces (right- and left-going electrons) by introducing the spinor 

We have neglected the small q-dependence of the matrix elements (i.e. the Fourier components of the impurity 
potential) in Eq. (|H). This is because we are dealing with quasi-one dimensional density wave systems whose Fermi 
surface consists of two almost parallel sheets at the points ±kp. Consequently, two relevant scattering amplitudes 
can be distinguished: the forward U(0) and the backward scattering parameter C/(Q), respectively. We believe that 
we can capture the essence of physics with this approximation. This is indeed the case, at least as long as only 
thermodynamics is concerned. Investigating the effect of pinning, and the sliding density wave however, requires the 
inclusion of the small q contribution . 28 ' 29 

Now, with all this, the electron thermal Green's function for the interacting system using Nambu's notation is 
defined as£ 

G a (k,iLu n ) = - dr(T T * a (k,r)*+(k,0)) H e^" T , (6) 
Jo 
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FIG. 1: Non-crossing self energy corrections due to impurity scattering. The solid line denotes the electron propagator G, 
while the dashed line is for the electron-impurity interaction. Dashed lines coming from the same cross represent successive 
scattering of the electron on the same impurity. 



where the Green's function is chosen to be diagonal in spin indices, and the subscript H refers to the total Hamiltonian 
H = Hq+Hi with which the expectation value should be taken. The propagator for the pure system without impurities 
is readily evaluated as 

nOfir ■ ^ _ iuJ n + $(k)/?3 + piRcA CT (k) - p 2 ImA tr (k) 

Lt (K, lUl n ) — 2 i t7T\2 i I A ft M2 ' 

+ £( k ) + l A <r( k )l 

where pi (i — 1, 2, 3) stand for the Pauli matrices acting on the space of left-, and right-going electrons. 

From now on we will drop the spin indices since they are irrelevant for our discussion on thermodynamic properties, 
and our results apply to both unconventional spin-, and charge-density waves as well. Furthermore, we choose the 
momentum dependence of the order parameter appearing in the zeroth order Green's function as A(k) = Ae 11 ^ sin(6fc y ), 
where the phase <j> is unrestricted due to incommensurability. Though the impurities we shall shortly introduce in the 
system are able to pin down the phase , 28 i 29 performing the usual impurity average restores translational invariance 
and consequently the phase degree of freedom. Thus, for brevity we fix the phase to be <j> — and without the loss of 
generality consider a real order parameter. As to the specific k-dependence, had we chosen a gap with cosine function 
or with the z component of the wave vector would not alter our results. 

Let us now turn our attention to the self-energy E(k, iu) n ) caused by the non-magnetic impurities. In diagrammatic 
language, the non-crossing approximation (NCA) retains only those diagrams where no two interaction lines cross 
each other. A simple subset of such diagrams are usually referred to as the ladder-type diagrams describing successive 
scattering processes on the same and single impurity atom. In case the particle lines appearing in these diagrams 
are not the bare G° but the dressed propagators G, then all the rainbow-type diagrams are included as well, and the 
NCA becomes full and self-consistent, see Fig. [T] The coupled equations valid for UDW are 3 ^ 

G- 1 (k,io;„) = Go 1 (k,^ n )-S, (8) 

y= U(0)-g(U(0) 2 -\U(Q)\ 2 ) fq s 
*l-2 g U(0) + g*(U(0)*- \ um 2y W 

where g = {2V)^ 1 TrG(k, iuj n ) is the local Green's function, rii is the impurity concentration, and according to the 
ignorance of small q-dependence of scattering amplitudes (see Eq. ((4]) and subsequent argument), the self energy turns 
out to be momentum independent E(k, iuj n ) = E(io; n ). We shall note here that the impurity concentration should be 
small enough in order for the NCA to be applicable and reliable. Quantitatively it means that the self energy should 
be much smaller than the bandwith W = 2vpkF, otherwise the contribution of crossing diagrams cannot be neglected 
and retaining the non-crossing class of diagrams only becomes inadequate. Moreover, at large enough concentration 
the standard independent average over the impurity positions on the lattice breaks down (for the impurity average 
performed on E, see Ref. [34]). We will see shortly however, that because of pair-breaking, we have a natural limit on 
rii which is of the order of 1%. Staying in this region of concentrations should be within the bounds. 
The ansatz-solution G of the closed pair of equations above reads as 



■ . t £(k)p 3 + A(k)pi 

G(1Wl) = ~ ^(k) 2 +A(k)3 ' ( 10 ) 



where, in contrast to conventional systems with constant gap, only the Matsubara frequency gets renormalizedi 34 ! 35 
The order parameter remains unaffected due to the unconventional nature of the gap with vanishing Fermi surface 
average. This is also known for quasi-one dimensional superconductors.-^ With the aid of the explicit expression in 
Eq. (fT0|) , g is calculated easily, and the renormalization condition for the Matsubara frequency follows as 

iuj n = iui n + Sfi — E, (11) 
9 = -No .f n k( f =], (12) 

V^n + a 2 Vy^ + A 2 J 
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where Nq = l/(%VFbc) is the normal state density of states per spin, K(z) is the complete elliptic integral of the first 
kind, and we have introduced a shift in the chemical potential Sfi by hand in Eq. pip . This we need in order to 
keep the particle number conserved in the system, since the presence of impurities on a macroscopic scale causes the 
particles to either flow in or out of the system depending on the sign and magnitude of the forward and backward 
scattering parameters. 

Impurity scattering from the theoretical point of view has been investigated extensively in the literature. These 
related works can be usually classified into one of two groups depending on the formalism used to capture the 
strength of the scattering process. However, both approaches are limiting cases of the formalism introduced here. In 
the second order Born approximation ) 9 ' 34 ' 35 ' 46 applicable for weak impurity scattering, the self energy simplifies to 
E = TiiU(0) — nig(U(0) 2 + \U(Q)\ 2 ), and S[i = rijZ7(0). On the other hand, for strong scattering the unitary limit is 
widely used ; 35 ' 47 ' 48 ' 49 where E = —rii/g and 8[i = 0. 

The aim of the present work is to make use of the full self energy defined in Eq. ©. Without making any 
constraints on the strengths of the scattering parameters U(0) and U(Q), we give a unified theory of the full non- 
crossing approximation. With this general approach we are able to cover the whole parameter space, to reobtain 
the well known formulas in the weak and strong scattering limits, and most importantly to provide exact results 
in the intermediate regime. As we have already mentioned, similar calculations based on a self-consistent T-matrix 
approximation have been performed for heavy fermion superconductors 3 - 9 -, and then later for d-density waves^ as 
well. In UDWs however this question is still open. 



III. THERMODYNAMICS OF IMPURE UDW 



Using the off-diagonal component of the Green's function in Eq. (flTJ|) . the gap equation, that determines the 
temperature dependence of the order parameter in the impure system, reads as 

A-^V Asin(6fc y ) 2 

N ^ Q 2 + £(k) 2 + A 2 sm(bk y ) 2 ' 1 ' 

where N denotes the number of sites. The effective coupling P favoring the sine dependence of A(k), and a factor of 
sin(6fcj,) in the numerator comes from the kernel of the electron-electron interaction.— Now we derive the equation 
that determines the transition temperature T c . Following Ref. we reformulate Eq. (fT5)) as 

(14) 




where E(z) is the complete elliptic integral of the second kind, and T c $ is the transition temperature of the pure 
UDW 4 s - We note that the basic Eqs. and (|1HI12|) determining the renormalized frequency as a function of the 
bare frequency provide unphysical solutions as well. The physically correct solution, that involves the right analytic 
properties of the Green's function, is chosen according to: sign(Re<D„) = sign(w„). In the region near the transition 
temperature T c , where A is small, we can expand the right hand side of Eq. (|14p . and we obtain the celebrated 
Abrikosov-Gor'kov (AG) formula 



-"{tJ =R '*{-2 + ^ +, ^fr) + *{-2)- (15) 

where ip(z) is the digamma function and 
Eq = RcE(A = 0) 

= 4n l f/(0)^- 1 (4 + N 2 tt 2 (U(0) 2 - |[/(Q)| 2 )), (16) 
E£ = ImE(A = 0) 

= -2 7 rsignK)n l iV ^- 1 [4(C/(0) 2 + \U(Q)\ 2 ) + n 2 N 2 (U(0) 2 - \U(Q)\ 2 )% (17) 
D = (4 + iV o V(t/(0) + |f/(Q)|) 2 )(4 + N 2 tt 2 (U(0) - |t/(Q)|) 2 ). (18) 

Up to this point the shift in chemical potential played the role of a constant parameter, the value of which is determined 
by particle conservation via the integrated density of states. However, the calculation and analysis of the one particle 
density of states is the subject of the next section, therefore we anticipate here the result 

SfM = jy . (19) 




Wo) 




Wo 1 17(0) 



FIG. 2: (Color online) Left panel: the shift in the chemical potential due to impurity scattering as a function of the scattering 
amplitudes. Right panel: contourplot of the scattering rate F/r c = 1 on the U(0) — |[/(Q)| plane for rii/n c = 1.15, 1.33, 2.08, 
2.56 and 4.76. For example in the case of rii — 2.56n c , the variation of the scattering parameters is restricted to the region 
inside the contour (blue) with the corresponding label. 



This result could have been presumed also from the AG formula in Eq. (JT5J), because in the Born and resonant limits 
the argument of ip is purely rea L 34 i 35 In Fig. [2] (left panel) we plot Sfi as a function of the scattering parameters. Note 
also, that the above formula is valid for normal metals with A = as well. With all this, we arrive to the final form 



1 2 2nT r 



(20) 



where we have introduced the scattering rate T = |Sq|. We see now, and shall point it out, that the AG formula is 
valid in general for arbitrary strength of the impurity scattering, and its applicability is not limited for the case of weak 
and strong scatterers only. The general form for T simplifies to the correct expression (Fi + T2)/2 in the Born limit, 
where Ti = nniNoU(0) 2 and T2 — irmNo\U(Q)\ 2 are the forward and backward scattering rates, respectively^ On 
the other hand, it gives back the result 2iii/(irNo) known for strong scatterers^ Its critical value, where T c vanishes, 
is given by r c = ttT c q/{2^) = -y/eAoo/4, with 7 = 1.781 being the Euler constant and Aoo stands for the pure zero 
temperature order parameter. One can easily see that T is bounded as a function of the forward and backward 
scattering parameters, but scales with n^. Consequently, this enables us to define a critical impurity concentration 
n c = it 2 N T c0 / (Aj) , below which (n, < n c ) the values of U(0) and |?7(Q)| can be arbitrary, thereby including the 
Born and unitary limits as well, but above which (rii > n c ) the amplitudes are constrained into a finite region on the 
U(0) — \U(Qi)\ plane around the origin. This is because otherwise the strength of the impurity scattering completely 
destroys the density wave condensate even at zero temperature. In the latter case apparently, only the Born limit is 
accessible. In Fig. [2] (right-panel) we show these regions and the corresponding confining contours of r/r c . Using the 
parameters of a-(BEDT-TTF) 2 KHg(SCN)4, namely T c — 10K, vp = 6 x 10 4 m/s^ and lattice constant in the chain 
direction a = 10 _9 m^ the critical concentration is estimated as n c — 0.01. Close to r c the transition temperature 
vanishes like^ 



27 Vr 



(21) 



Next, we proceed with the analysis of the zero temperature gap maximum. For brevity, in this paragraph A 
A(T = 0). Similar calculations that lead to Eq. ([I4])jlui give in our case 

ro ° dto 



U) 



(22) 



where uj and il> are, respectively, the bare and the renormalized Matsubara frequencies (we omitted the subscript n), 
and 
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The expression on the right hand side of Eq. (|22|) , though is not adapted well for the numerical solution, is however 
convenient for further analytical calculations. Namely, it enables us to determine the behavior of the order parameter 
close to the critical scattering rate T c . As the transition temperature and the gap maximum are closely related 
quantities, we expect similar tendency found for T c in Eq. (|2ip . Indeed, we obtain 



4Te 

V3 



2ttRcA 

~9~T7 



-1/2 



1U3|ES 



(24) 



where A = NodT,(go)/dg, and go = g(A = 0,u> n > 0). As we have just pointed out, for numerical purposes Eq. ([22 
shall be reformulated, and the following form is obtained 



2 V Aoo 



Co 



dzF{z) 



2 
A 



dz F(z)lm 



<9£ 



2C> 



d\Im<F(z) 1 



i <9E 
A ~dz~ 



(25) 



where z = u>/ A, Co = Cq + iC ' is the decomposition of Co into real and imaginary parts, and in the last integral 
z = C + iACq . The Co parameter is the value of <D„/A at zero frequency as usuah 34 i 35 ' 47 Note, that due to the true 
complex nature of the self energy X with finite real part, this generalized Co acquires a finite imaginary part leading 
to the appearance of the third term on the right hand side of Eq. (f25|) . We present the result of the numerical solution 
of Eq. (|2"5"1) in Fig. [3] (left panel) for fixed impurity concentration iii/n c = 2.56, that exceeds the critical threshold n c . 
We shall emphasize, that the impurity concentration and the scattering parameters enter in the expression for A not 
through r, as is the case for (U)DWs^ ' 34 ' 35 and superconductors^ in the Born and unitary limits and also for T c in 
the AG formula (see Eq. (|20|) ). but enter through £ and Sp. In this respect, the plot shows the dependence of the 
zero temperature order parameter on the scattering amplitudes above the U(0) — |f(Q)| plane, with m fixed. One 
can easily see, that the contour where A vanishes and T takes on the critical value r c (see Eq. (|24I0 . is just the one 
appearing in the right panel of Fig. [21 with the corresponding label 2.56. 

The dependence of T c on the scattering amplitudes is qualitatively the same as that of A(0) in Fig. 02 therefore is 
not presented here. However, the A(0)/T c ratio might be of interest, and is shown in the right panel of Fig. [3] It 
is clear from the figure, that the ratio increases with increasing strength of impurity scattering. The lower bound is 
obviously the pure value 2n/( , yy/e) — 2. 14. 43 Moreover, the maximum is reached on the confining contour 



max 



A(0) 



8 Tie 

5 rii 



-1/2 



(26) 



Here we have used Eqs. ([2"T]) and |24|) . Now it is clear that as n* increases, the overall maximum of the gap to T c ratio 
is y/8/5ir = 3.97. This result has been obtained previously for quasi-one dimensional anisotropic superconductors 
based on a second order Born calculation^ 5 - 

Close to T c , A(T) vanishes in a square-root manner as does usually in mean-field treatments. In the series expansion 
of Eq. (TT4"]) for small A, the term of order A 2 provides 



8(27rT c 



p) 



1 RcA 
12 T c 



T 



(27) 



where ip^(z) is the nth polygamma function and p — T/(2nT c ). The term (12T c ) _1 ReA in the denominator simplifies 
to p/3 in the Born^ and to — p/3 in the resonant unitary 4 ^ limits. 

Now we derive an expression for the order parameter close to absolute zero. The alternative formulas for the gap 
equation in Eq. (|13p presented so far are not really convenient for our present purposes, thus with the usual contour 
integration technique we rewrite it as 



4 V Aoo 



du> ( to _ 
— <^ —Re 
A I A 



oj J \ to 



- Re 



u 



Rc 



in* 



U 



(28) 



where A means A(T), and u is the analytically continued value of i£j n , namely u = i£j n (iuj n — > to + iO). In addition 
f(u>) is the Fermi function and the subscripts ±w mean that u has to be evaluated at the corresponding frequency. 
This form will also be used for the evaluation of the low temperature specific heat shortly. We note that in the pure 
case, Eq. (|28[) simplifies considerably, and one can recognize the result found for pure d-wave superconductors. 51 This 
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FIG. 3: (Color online) Left panel: the zero temperature order parameter A(0) as a function of the scattering amplitudes for 
fixed impurity concentration rii/n c = 2.56. Right panel: the A(0)/T c ratio versus the scattering amplitudes for the same 
concentration as in the left panel. It is clear that the presence of impurities enhances this ratio from the pure result 2.14. For 
both panels, the contour of the domain on the U(0) — \U (Q)| plane is the blue curve in Fig. [2] right panel. 



is not surprising at all, as the similarity between the thermodynamics of UDW and d-SC are well-knowni 34 i 35 i 43 i 47 
From Eq. (|28|) we finally obtain 



where the dimensionless parameters a and b are 



--*H(^)">(^y 

4 f° dcu [ , r J dYX~ x d , /A(0) 
- / ^tRc US - 5u 1 + — —E' — ^ 



(30) 



h =~ -T-7-rRc < (S — <5/i) 1 + — I -E' • (31) 



7T 



A(0) y ^ \ du) du 



In Eqs. (|30j|31[) uq = 2,A(0)Co, and the prime means differentiation with respect to the argument. With cumbersome 
calculations one can show that Eq. (|2T))) gives back the known expression in the weak scattering Born limit i^ 4 - The 
most important message of the above formula is obviously the fact that the presence of impurities changes the pure 
T 3 power law behavior— and transforms it to a faster T 2 decrease. 

We shall close this section with the presentation of some relevant thermodynamic quantities, such as the grand 
canonical potential f2, the entropy S and the specific heat Cy . Our starting point is the well-known formula for 
due to Pauli, that is 

f 1 dX 

n-n Q = —(\H int ), (32) 

where Hi nt is the interaction causing the phase transition. This formula gives us the thermodynamic potential 
difference <5fi between the normal and the DW phase. After lengthy calculations we obtain 











'—) 


+ 2 / duiuj( 




,Aoo/ 


J — oo \ 



where A means A(T) everywhere. Moreover N(ui) and N p (ui) are, respectively, the quasiparticle density of states 
(DOS) in the impure and pure^ 3 - systems, and finally N s (uj) is the symmetrized DOS: N s (uj) — (N(u) + N(—lu))/2. 
We will come back to the analysis of DOS in the next section, and we will see that it is not an even function indeed. 
We shall emphasize that the order parameter which is hidden in N p (lo) and determines the energy scale, is that of 
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FIG. 4: Left panel: the residual density of states at the Fermi level for fixed impurity concentration rii/[NoA] = 0.3. Right 
panel: the relative jump in the specific heat at the transition temperature T c normalized to the pure value as a function of the 
scattering amplitudes for fixed impurity concentration ni/n c = 2.56. 



the impure system. We also note that Eq. (J33J) is valid for all temperatures. In obtaining the above formula, we 
repeatedly made use of Eq. (|2"B"|). and the following important relation 

Now, one can easily check, that the temperature dependence through A does not contribute to the entropy since the 
gap equation ensures dSfl/dA = 0. Consequently, in an UDW we simply get 

/>oo 

S = 4 / dujN s {uj) (ln(l + e-P") + , (35) 



o 

which is the standard expression of a normal metal, except the DOS is that of the impure unconventional density 
wave. Note that due to the asymmetric feature of the DOS with respect to the Fermi energy, here N s (lu) shows up. 
In those systems however, where the particle-hole symmetry is not violated, this feature remains hidden obviously. 49 
From this, the low temperature specific heat is Cy = (2ir 2 /3)N(0)T. Here N(Q) stands for the residual density of 
states on the Fermi level and it is shown in the left panel of Fig. 0J 

Around T c the jump in the specific heat, normalized to the pure result, follows as 

21cf3) Jc (i-^(i + P)) 2 m) 

AC V0 ^ W T c o|^(I+p) + -L^'"(i+p)' W 

where ACyo = 167T 2 N T c0 / (21( (3)), and it is shown in the right panel of Fig. 2] for fixed impurity concentration 
rii — 2.56n c . This result, just like Eq. ([!??)) , also becomes simpler in the case of weak and strong scattering limits, 
(12T c ) _1 ReA — > ±p/3, and we reobtain the familiar results^ Note that in Ref. |H a factor of T c /T c q is missing, but 
is corrected by the same authors in a separate paper.— 

IV. DENSITY OF STATES IN UDW 

With the knowledge of the Green's function (see Eq. (jTUJ)) , the quasiparticle density of states can be calculated as 

k 

= — iVo Im | K (—==] ) = -No ReK ( — ) , (37) 
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FIG. 5: Left panel: the quasiparticle density of states as a function of reduced energy for different values of the scattering 
amplitudes at fixed impurity concentration n t /[iVoA] = 0.4. U(0)/\U(Q)\ = 2 and N \U(Q)\ = 0, 0.2, 0.4, 0.6, 1, 1.5, and 20 
for curves with increasing N(0). Right panel: the density of states enlarged around the Fermi energy for m/[NoA] = 0.02, the 
ratio of the two scattering parameters is the same as in the left panel and iVo|t/(Q)| = 0, 0.6, 1, 1.5, 3, 5, 10, and 50 for curves 
with increasing N(0). 



where u = iuj n (iuj n — > u> + iO) is the analytically continued value of the dressed Matsubara frequency. It is determined 
self-consistently from 



u = lo + Sfi — 



(38) 



and as we have already pointed out before, the physically correct solution is the one that ensures the right analytic 
properties of the retarded Green's function, that is Imu > 0. 

Now, we turn our attention to the shift in the chemical potential appearing first in Eq. and here in the 

analytically continued version Eq. (|3"8|) as well. We have noted in the previous section also, that S/i was introduced 
by hand in order to maintain particle conservation. Its actual value is to be determined from the integrated DOS as 



(1uj{N(lu) ~ N ) = 0. 



(39) 



The second term in the integrand is the constant density of states of the pure metal, but it could equally be the 
DOS of the pure UDW. This is because the metal to UDW transition without impurities conserves the particle-hole 
symmetry and the even parity of the DOS, and as such, conserves particle number. The numerical solution of Eq. (|39p 
confirms our expectation (based on the AG formula, see Eq. lfl"5])). and we indeed find = Eg. 

With all this, we are now able to numerically evaluate the quasiparticle DOS defined in Eq. ([3T)) . and the results 
are shown in Fig. [5] for positive, i.e. repulsive forward scattering. Apparently, the main feature of these curves is the 
violation of particle-hole symmetry, similar to that found in d-density waves 4 ^. Namely, strictly speaking they are 
not even functions of energy. However, in the weak and strong scattering limits we reobtain even parity ] 34 ' 35 ' 45 In 
addition, it is clear from the figures, that the logarithmic singularities at u> = ±A characteristic to UDWs in general 
are smoothed out due to the presence of scatterers. In the case of repulsive (attractive) scatterers the upper (lower) 
peak is larger. An other well-known effect of impurities is the finite residual DOS at the Fermi energy, N(0) > 0. For 
small forward and backward scatterings it disappears exponentially, it is shown separately in the left panel of Fig. [4] 
The right panel of Fig. [5] demonstrates the behavior of DOS on a smaller energy scale around the Fermi energy. For 
small rii and large scattering amplitudes, we come up with the familiar result of the resonant limit ) 35 ' 47 ! 53 that the 
density of states exhibits an island of localized states around the Fermi energy superimposed on the usual gapless 
density of states of the pure system, which is manifested in the nonmonotonic nature of the DOS close to w = 0. 

After some algebra, the low-energy behavior turns out to be linear, and it reads 



N(w) = JV(0) + -7V Im 




(40) 
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where the argument of the elliptic integrals is This is to be contrasted with the lo 2 contribution of the 

extreme limits. On the other hand, at high energies the density of states reaches the normal state value as 




This is what we find in the second order Born calculation,— except that now Y is defined by the general formula 
T = | Eg | (see also Eqs. f|lTtil8[) V It is valid for arbitrary strength of impurity scattering parameters. 



We have studied the effect of non-magnetic impurities in quasi-one dimensional unconventional charge and spin 
density waves within mean-field theory. The focus has been on the thermodynamic response of the system. As the 
impurities are of the ordinary non-magnetic type, the spin degrees of freedom in the condensate play no role, leading 
to the same conclusions for both types of UDW. As opposed to the usual treatment of impurity scattering in density 
waves and superconductors (i.e. the Born and unitary limits), we do not put any constraint on the strength of forward 
and backward scattering amplitudes, but make use of the full self energy obtained in the self-consistent non-crossing 
approximation (NCA). Namely, the applied proper self energy includes all ladder-, and rainbow- type diagrams to 
infinite order. In this regard, our NCA is full. We than repeatedly took advantage of the fact, that our generalized 
formulas contain the familiar results of the weak and strong scattering limits, and we reobtained them indeed. On 
the other hand, our results provide explicit expressions in the intermediate regime as well. 

We started our investigation with the transition temperature and the order parameter affected by the scatterers. 
We have found that, with a suitable redefinition of the scattering rate compared to either the Born or the unitary 
expressions, the Abrikosov-Gor'kov formula holds in this general case as well. Its validity can be therefore extended 
to arbitrary scattering amplitudes. The analysis of the AG formula allows one to introduce a critical impurity 
concentration, below which (rii < n c ) the impurity potential, whatever strong it is, cannot break down the density 
wave condensate. However, if rii exceeds the critical threshold n c , the order of which is estimated around 1%, the 
scattering amplitudes cannot be arbitrary large and only the Born limit is accessible. An interesting consequence 
of impurity presence is the appearance of a shift in the chemical potential 5/i, the inclusion of which is necessary in 
order to conserve particle number. We gave an explicit and exact expression for Sfj, as a function of the scattering 
amplitudes and concentration, that is valid in a normal metal as well. The overall maximum of the zero temperature 
gap to T c ratio was found to be the same (-\/8/57r = 3.97) as in anisotropic quasi-one dimensional superconductors. 




The most eye-catching feature of impurities is definitely the violation of particle-hole symmetry manifested in the 
asymmetric quasiparticle density of states. Strictly speaking, the DOS is not an even function of energy with respect 
to the proper Fermi energy, that is the pure Fermi energy shifted by 5(1. The asymmetry is most apparent around 
lu = 0, where the low frequency behavior is linear as opposed to the quadratic form of the extreme limits. This also 
involves that the minimum of DOS is not at the Fermi energy. In addition, the peaks around lo = ±A, the residues of 
logarithmic singularities of the pure DW, become asymmetric as well: in case of repulsive (attractive) scatterers the 
upper (lower) peak is larger. In unconventional DWs, at any finite scattering strength the valley of the DOS around 
the Fermi energy is filled in, leading to normal electronlike behavior very close to absolute zero, but the reduced 
density of states compared to the normal state signals the effect of the condensate. Taking the weak and strong 
scattering limits of the DOS, one regains particle-hole symmetry, and can easily reobtain the familiar lincshapes of 
the Born and unitary limits, in particular for instance the island of localized states around oj = in the latter case. 
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